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Abstract 

Treatises on vibrations devote large space to study the dynamical behavior of an 
elastic system subject to known external tractions. In fact, usually a "system" is 
not an isolated body but it is part of a chain of mechanisms which disturb the 
"system" for example due to the periodic rotation of shafts. This kind of problem 
has been rarely studied in control theory. In the specific case we shall study, the 
case of a viscoelastic string, the effect of such external action is on the horizontal 
component of the traction, and so it affects the coefficients of the corresponding 
wave type equation, which will be time dependent. The usual methods used in 
controllability are not naturally adapted to this case. For example at first sight it 
might seem that moment methods can only be used in case of coefficients which 
are constant in time. Instead, we shall see that moment methods can be extended 
to study controllability of a viscoelastic string subject to external traction and 
in particular we shall study a controllability problem which is encountered in the 
solution of the inverse problem consisting in the identification of a distributed 
disturbance source. 
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1 Introduction 

Elastic and viscoelastic systems have been widely studied from the point of 
view of controllability (see for example [9j[I7J|20] and references therein), the 
simpler "canonical" case being the string equation 



The coefficients c(x) = l/p(x) depends on the physical properties of the 
string i.e. the density, and P is the orizontal traction in the string which, 
for small oscillations and isolated systems, is constant in time and space 
(while c may depend on x since the density may not be constant.) It is 
usually explicitly assumed in control problems that, a part the action of 
external controls which does not change the horizontal component of the 
stress, the string is "isolated" i.e. that P does not depend on time. Control 
problems for the string equations under external traction have been rarely 
studied, see [TJ [2j [3]. In fact, in general the system is part of a chain of 
mechanisms which produce "perturbations" which influence the horizontal 
traction in the string, often in a known way, for example due to the rotation 
of shafts. For this reason, books on vibrations devote large spaces to the 
study of the behaviour of elastic systems under known external tractions, see 
for example [TJ [15] . 

In this paper we are going to study a viscoelastic string subject to exterior 
known tractions. The string is also subject to a known signal a(t) which 
enters trough an input operator b which is unknown. For example b may 
be the characteristic function of a certain set and our final goal will be the 
identification of the input operator b using boundary observations. So, the 
problem we are going to study boils down to the following equation 



where P(t) represents the effect of the variable external traction and w = 
u>(£, t) with t > and £ G (0, it) (independent variables in w are not indicated 
unless needed for clarity). 





Remark 1: Writing the equation in the form (QP implies that the string is at 
rest for negative times. If not, a known additional terms will appear in the 
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right hand side of the equation, which has no influence for the arguments 
below. So, we put it equal to zero, i 

Initial and boundary conditions are homogeneous (of Dirichlet type) 

w(£,0) = 0, w t {£,0)=0, w(0,t) = w(n,t) = 0. (2) 
The observation is the traction at x — 0, i.e. 

rt 



-P(t) 



c(O^(0,t)+ f M(t-s)c(£)wt(0 t s) ds 
Jo 



We shall assume that P(t), M(t) and c(x) are smooth and thatP(t) and 
c(x) are strictly positive (see below for the assumptions) so that from this 
"real" observation we can reconstruct the function 

V (t)=wz(0,t) (3) 

which will be used in the reconstruction algorithm. 

In the purely elastic case, this kind of inverse problem has been studied 
using control ideas in [18] (see also [6j [14] ) where it is proved that this inverse 
problem (in the purely elastic case) depends on the solution of a boundary 
control problem. So, it is easily guessed that the control problem has to be 
studied first also in the case that the traction varies in time. This is the 
second main subject of this paper which is studied first. 

The problem we need to study first is the controllability of the deformation 
w(-,T), at a certain time T, when b = 0, the initial conditions are zero but 
the boundary conditions are 

f(t) 

w M = 4$(fi> w(n,t)=0. (4) 

Here f(t) is a square integrable control (the denominator is introduced solely 
for convenience). So, the control problem we shall study first, and which 
accupies the most part of this paper, is as follows: given any prescibed target 
W(£) G L 2 (0,7r), we must prove the existence of a square integrable control 
f(t) G L 2 (0,T) such that w(£,T) = W(£). Furthermore, we shall see that T 
does not depend on W. 

The organization of the paper is as follows: the control problem is studied 
in section |3] while the identification problem is studied in section HI Assump- 
tions and preliminaries are in Section [2j 
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1.1 Further references 

Control problems for viscoelastic materials have already been studied under 
several different assumptions. We cite in particular [H [5] which seems to con- 
tain the most general results. The paper [8] studied the case M = M(t, s) but 
it is explicitly stated that the arguments require constant density. The pa- 
per [5] studied the case M = M(t, s, x) when the control is distributed (from 
which controllability results for boundary control systems with the control 
acting on the full boundary can be derived). These papers study also the 
multidimensional case, but the results does not seems to be "constructive" 
since they depend on compactness arguments (paper [5], which explicitly re- 
quires constant density) or Carleman estimates (paper [S]). In this respect 
we cite also [IT], where controllability under boundary control is proved for 
multidemensional systems. The arguments are presented in the convolution 
case and constant density, but are easily extended to non constant density. 
The proof in this paper being based on compactness arguments is not con- 
structive either. The paper j3], being based on an extension of D'Alambert 
formula (see last section) is more constructive in spirit, but the elastic oper- 
ator is assumed to be time invariant. The papers [121 H3J HHj instead study 
the control problem for a viscoelastic string, using moment methods, so that 
a representation formula for the control can be given (see in particular [121 
Section 4] and [101 Section 4]). Our goal here is the extension of these re- 
sults to a viscoelastic string subject to nonconstant traction, and to show 
that the results can be used to solve a source identification problem. Source 
identification problems using control methods have been solved for elastic 
materials with constant traction first in [18] (see also [6]) and the method 
has been extended to materials with memory in [13]. When the traction is 
not constant we need a different idea, presented in section HI 

2 Assumptions and preliminaries 

The assumptions in this paper are: 

1. P(t) is continuous and strictly positive, P(t) > po > for every t > 
while c(£) G C 1 (0,7r) is also strictly positive: c(£) > c > for every 
£e [0,4 
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2. The kernel M(t) G Wj c 2 (0, +00). Hence its primitive N(t) 

N(t) = 1 + / M(s) ds 



has three derivatives (N"'(t) is locally square integrable) and N(0) = 1. 

We introduce the following selfadjoint operator A: 

dom A = H 2 (0, tt) n ^(0, tt) , A0 = (c(0^(0){ • (5) 

It is well known that this operator has compact resolvent, and that it has a se- 
quence of eigenvalues {— A^} n >i, with the following asymptotic estimate [HI 
p. 173)] 

\ n = n + — , |#„|<M (6) 
n 

A corresponding sequence of real normalized eigenf unctions <fi n {V) has the 
following estimates, see [161 P- 176-177]: 

M n (0 



0n(O = \/~ sin< if(£)cosn£ 



7r v n " I n 2 



_[sinr* + ^ (l^(OI<M, \K'(0\<M, 
A n /' IIU0I<M, \M n (0\<M, [() 

#,(0 = yf {ncosne + fcn(0} > |fcn(C)l<M. (8) 



The previous formulas have been taken from [16] but this reference has 
n + 1, with n > 0, where we wrote n (and we intend n > 1). Sharper 
estimates are known, see [161 formula (82) p. 173 and (85) p. 177] but we 
don't need this. 



2.1 Riesz systems 

The arguments of this paper are based on the theory of Riesz systems. We 
recall the definition: a sequence {h n } in a separable Hilbert space H is a 
Riesz basis when it is the image of an orthonormal basis of H under a linear 
bounded and boundedly invertible transformation. The key property is that 
every h e H ha a unique representation with respect to the Riesz basis {h n }, 



h=) j a n h n 
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and {a n } G I 2 . In fact, there are positive numbers m and M such that 

m\\h\\ 2 < ^2\a n \ 2 < M\\h\\ 2 . 

If a sequence is a Riesz basis of its closed span (possibly not equal to H) it 
is called a Riesz sequence. 

A property we shall need is that if {h n } is a Riesz sequence then the series 
Y^Oi n h n converges if and only if {a n } G I 2 , and the convergence is uncondi- 
tional. Furthermore we need two perturbation results from [T9~| Ch. 1]. The 
first is a Paley- Wiener theorem, adapted to Hilbert spaces and orthonormal 
bases. It states that if 

^2\\h n -e n \\ 2 < 1 

({e n ,} orthonormal) then {h n } is a Riesz basis. In fact we shall use the 
following corollary: 

Corollary 2: Let {e n } be a Riesz sequence and let the sequence {h n } satisfy 

W h n — e n || 2 < +oo (9) 

then there exists a number N such that {h n } n> N is a Riesz sequence too. 
This in particular implies that 

/] Ot n h n 

converges in the norm of H if and only if {a n } G I 2 . 

We stress that the sequence {e n } in Corollary [2] need not be orthonormal. 
Condition Q does not imply that {h n } is a Riesz sequence but 

Theorem 3 (Bari Theorem): If the sequence {h n } satisfies condition and 
if furthermore 

J2a n h n = {a n } = (10) 

then {h n } is a Riesz sequence. 

The additional condition (fTUj) is called w-independence. Note that the 
convergence of the series in (fTUj) is in H so that it must be {«„} G I 2 because 
{h n }n>N is a Riesz sequence. 
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2.2 Representation of the solution 

Now we consider Eq. ([!]) with zero initial conditions but both the affine term 
and the boundary condition (j3J). We need a definition/representation formula 
for the solutions. We confine ourselves to a formula for w(x, t) and we ignore 
velocity, since only w(-,T) is needed for the identification problem. So, we 
prefer to write Eq. (DQ) as a first order equation, integrating both the sides. 
We get 



w t 



N(t - s)P(s) (c(0^) £ (s) ds + b(£)g(t) , g(t) 



a(s) dS. 
(11) 

We multiply both sides of (ITT]) with 4> n (£) and we integrate on [0, w]. We get 

<(t) = -A* f N(t - s)P(s)w n (s) ds + #,(0) f N(t - s)f(s) ds + b n g(t) 
Jo Jo 

(12) 



where 



w n (t) = / w(£,t)<f> n (£) dx, b n = (b,(f) n ) = 
Jo 

Let us fix any time T > 0. We want a representation formula for w n (T). 

Let 

Q(t)=P(T-t) 

and let z n (t; T) solve 

z' n (t; T) = -\ 2 n Q(t) [ N(t - s)z n (s; T) ds , z n (0) = 1 . 
It is simple to check the following representation formula for w n (T): 



w n {T)= [ T f(T-s) k(0) / 

^0 ^0 



0'„(O) / N{s-r)z n {r-T)dr 



ds 



+b n / z n (T - s;T)g(s) ds . 
Jo 

so that we expect the following definition of the solution w(£, T): 

+oo „t 

w{£,T) = / f(T-s) 

n=l Jo 



<t>' n (0) / 7V( s -r)^(r;T) dr 



+oo „i 

+ y24>n(£)K z n (T - s;T)g(s) ds . 

n=l J ° 



(13) 



ds 



(14) 
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The following Lemma justifies formula (I14p : 

Lemma 4: The series in p4\ ), as functions ofT, converge in C([0, K}; L 2 (0, tt)) 
for every K > 0. 

The proof of this Lemma requires a preliminary study of the functions z n (t, T) 
and can be found in the appendix. 

3 The control problem 

In this section we study the control problem needed to identify the source 
term b; i.e. we consider Eq. ([I]) with b = and initial boundary conditions 

f(t) 

w(£,0) = 0, w t (£,0) = 0, w(0,t)= / n }J, A , w(n,t)=0. (15) 

c(0)P(t) 

As we already said, we want to identify a time T such that for every 
W(£) G L 2 (0,7r) there exists a control f(t) such that the corresponding 
solution w(t) satisfies w(T) = w(£,T) = W(£). 

It is convenient to perform several transformations first. In particular, as 
we are interested to the controllability solely of the deformation w(-;T) and 
not of the velocity, we prefer to rewrite Eq. (Q]) in the form of a first order 
equation, but we note explicitly that the technique we use here can also be 
used to study the controllability of the pair (w(-,T),w t (-,T)) as done in [10] 
in the case that the external traction P and the density are constant. 

This section on controllability consists of three subsections: first we re- 
duce the control problem to a moment problem. Then we present the trans- 
formations which will simplify the problem uder study, in subsection 13.11 
Then we prove that a certain sequence of functions is a Riesz sequence in 
L 2 (0,S) wher S is a suitable number to be identified (in Theorem [TO]) and 
this will prove controllability, as seen in the subsection 13.11 The final result 
of controllability is as follows (compare ( 1321) ): 

Theorem 5: Let T solve 

3 ^/Q{s) ds = [ ° ^P{T-s) ds = ir 
Jo 

and let T > T . Then, for every £ G L 2 (0, tt) there exists a boundary control 
f G L 2 (0,T) such that the corresponding solution w(t) satisfies w(T) = £. 
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Note that the number T exists since 

/ y/P(T-s) ds>Tp . 
Jo 

The computation of w(t) for every t makes sense tanks to Lemma HI 

Remark 6: In this Section the final time T is fixed once and for all. So, 
instead of the notation z n (t,T) we shall use the simpler notation z n (t), i.e. 
we put z n (t) = z n (t,T) (T fixed). We shall return to the complete notation 
in Section [^] and in the Appendix, i 

3.1 Reduction to a moment problem 

We fix a certain time T and we consider the solution w(£,T), at this fixed 
time, when b = 0. We have, from (Tl3|) and (fT4|) : 



w 



(£,T) = X>(0 ^ f(T-s) k(0) / 

n=l ^ L ^ 



<j)' n (0) / N(s-r)z n (r)dr 



ds. (16) 



As we noted, z n {t) instead of z n (t; T) since T is fixed. 

Equality ffTBT) shows that w(C,,T) = W(£,) if we can solve the moment 
problem 



f(T~s) 



0n(°) / N(s-r)z n (r)dr 



ds = W„ 



(17) 



where 



Note that the transformation of L 2 (0,n) 3 W t— > {W n } is an isometric 
isomorphism of L 2 (0,ir) and I 2 . 

Our main result is that this moment problem is solvable and that the 
solution / of (|17p (which has minimal L 2 -norm) depends continuously on £, 
since we have: 
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Theorem 7: Let T = T be as in Theorem^ and T > T . 
The sequence 

{<f>n(0)J o N(t-r)z n (r)dr} (18) 

is a Riesz sequence in L 2 (0,T). 

In order to prove the theorem, we need to perform several transforma- 
tions. 

3.1.1 Transformations 

Again we recall that we work with a fixed time T and < t < T so that 
z n (t; T) is simply denoted z n {t) but it is important for the Appendix that we 
don't forget the dependence on T . 

In this section, we transform the equation of z n (t). 

First we introduce 

H(t) = -N'(0)t-logQ(t), C n ( t ) = e H ^z n (t) = e- N '^ t ^ r) z n (t). (19) 

We shall see below the reason for this misterious looking transformation 
which does depend on T since Q(t) = P(T — t). 
The function ( n (t) verifies 

C n (t) - H'(t)( n (t) = -\ 2 n e H ^Q(t) f N{t - s)e- H ^Us) ds 

Jo 

and so also 

C - H'{t)C + (X 2 n Q(t) - H"{t)) Cn = G(t) (20) 

where 

G(t) = -X 2 n [H'(t)Q(t)e H ^ + e H ®Q'(t)] f N(t - s)e~ H ^Us) ds 

Jo 

-\ 2 n e H ^Q(t) f N'(t - s)e- H ^Us) ds . (21) 
Jo 

Now we apply Liouville transformation to Eq. fl20p . which requires the 
introduction of two auxiliary functions, which will be determined later on 
(see [HH p. 163]): we introduce the function ( n (t) given by 

a(t)Ut) = Ut) (22) 
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(here a(t) is the first auxiliary function). Then we have 
a(t)C(t) + [2a'(t)-H'(t)a(t)}C n (t) 

+ [{X 2 n Q(t) - H"{t)) a(t) - H'(t)a'(t) + a"(t)] = G W ■ 

The second trasformation is a transformation of the time variable: we intro- 
duce 

x = L(t) 

where L(t) is still unspecified (but we shall see that we can choose L(t) 
strictly increasing and its inverse transformation will be denoted M(t).) We 
introduce Y n {x) defined by 

( n (t) = Y n (L(t)) . 
Then we have the following equation for Y n (x): 

a{t) \L\t)f Y"{L{t)) + [a{t)L"{t) + (2a'(t) - H'(t)a{t)) L\t)} Y'(L(t)) 
+ [X 2 n Q(t)a(t) + (a"(t) - H"(t)a(t) - H'{t)a\t))] Y(L(t)) = G(t) . 

We must be precise on this point: here and below notations like Y'(L(t)) will 
denote the derivative of Y(x) computed for x = L(t). 

Now we relate the functions L(t) and a(t) so to have the coefficient of 
Y\L{t)) equal to zero, i.e. we impose 



(L'(t))' = -L'(t) 



a[t) 



(which is legitimate since we shall choose a(t) ^ 0.) In order to satisfy this 
condition we choose 



L'{t) 



e H(t) 



a\t) ■ 

Note that L(t) turns out to be strictly increasing. Furthermore, we choose 



* e H(s) 



so that we have also L(0) = and L(t) > for t > 0. 
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Once L(t) has been chosen as above, we see that the equation of Y n (x) is 
1 



fl 3(t) 

= G(t) . 
Finally we choose 



e 2H(t) Y:m) + [x 2 Q{t) + (aft{t) _ H „ {t)a{t) _ H / (t)a / (t)) ] Yn (t) 



a(t) = e 11 ^ 2 



(note: it is strictly positive) and we get the final form of the equation for 
Y n (t): 



Y"(L(t))+ K + V(t) Y(L(t)) 



e -H(t)/2 

WF 1 



G(t) 



(23) 



where 



e -H(t)/2 

V(t) = —. ,„,. [a"{t) - H"(t)a{t) - H'(t)a'{t)} . 



Q{t)V* 
The definition of a(t) gives 



L\t) = #j 



(24) 



Let M(x) be the inverse function of L(t), defined on the interval [0, L(T)] = 
[0, S] and 

V(x) = V(M(x)) . 

Using this and the explicit form of G(t) in ( 12T1) we see that Y n (x) satisfies 
the following integro-differential equation on the interval 



< x < S = L(T) : 



Y"(x) + {Xl + V(x)) Y n {x) 

= -\i e mm)/2 [h'{m(x))q x i\m{x)) 

i-M{x) 

+Q'(M(x))Q- 3/4 (M(x))] / N(M(x) - s)e"^ (s) C„(s) ds 

Jo 

pM(x) 

-X 2 n e Hmx))/2 Q 1/4 (M(x)) / N'(M(x) - s)e- H{s) ( n (s) ds 
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We make the substitution s = M(r) (we recall that M(r) is the inverse 
function of L(s)), we use < r < x when < s < Mix) and 

— — M(r) = - , . = — , where s = Mir). 

dr 1 ; L'{s) ^/Q(s) 1 ' 

Then, we have the following equation for Y n = Y n (x): 

Y: + V(x)Y n + X 2 n Y n = -X 2 n / A(x, r)Y n {r) dr (25) 

Jo 

where 

A(x,r) = e H(M(x))/2 Q 1/4 (M(x)) {[H'{M{x)) 
+Q'(M(x))Q- 1 (M(x))] N(M(x) - M(r)) 
+N'(M(x) - M(r))} e" H(M(r))/2 g- 3/4 (M(r)) . (26) 

The values Y n (0) = y and Y^(0) = y\ are easily computed: 

y = Y n (0) = e H ^^/Q{0) = Q(0r 1/4 = P{T)^\ 
Vi = K(0) = -\Q(0y 7/i {2N'{0)Q(0) + Q'(0)} 

= ^P(T)- 7 ^{P'(T) -2N'(0)P(T)} . 

In the following, we don't need these explicit expressions, but we shall 
need the following facts: 

• the initial conditions yo and y\ do not depend on n; 

• the initial condition y is strictly positive (we shall need that it is 
different from zero.) 

• the initial conditions y and y n do depend on T, equivalently on S = 
L(T). 

In the appendix we shall use the fact that when T e [0, K] then there 
exists a number M (which depends on K) such that |f/g(T)| < M, 
\yi(T)\ < M for every T G [0, K]. 



Viscoelastic string under external traction 



14 



We introduce the notation 



g n (x) = y cos \ n x + —y x sin \ n x . 



(27) 



Using (|27|) and ( 125|) we can write the following Volterra integral equation 



for Y n (x): 



1 f x 

Y n (x) = 9n{x) - — / sin A n (x - s)V(s)Y n (s) ds 

*n JO 

px PS 

— A n / sinA n (x — s) / A(s, u)Y n (u) du ds 
Jo Jo 



= 9n(x) - < / A(x,s) + —sin X n (x - s)V(s) 
Uo L a « 

+ / cos A„(x — s)A(s, s)F n (s) ds 
Jo 

PX PS 

+ / cosA„(x — s) / A i(s, i/)Y^(i/) dv ds 



Kfs) ds 



In this and following formulas, we use the comma notation for the deriva- 
tive. Hence, A^(x, s) is the derivative of A(x, s) respect to the first variable. 
When x = r the brace in fl26l) is 



{H'{M{r)) + Q\M(r))Q~\M(r)) + N'(0)} = 

thanks to the choice ffT^l) for H(t). The reason for the choice of the expo- 
nential e H( -^ is precisely this: 

A{s,s) = 0. 

So, elaborating further the integral equation for Y n (x) we have 



Y n (x) = g n (x 
1 

JO 
1 rX 



AJx, s)YJs) ds 



+— / sin X n (x - s) 



[ A A1 (s,v)Y n (u) dv 
Jo 



ds 



AJx,s)Y n (s) ds 



An Jo 



YJu) 



A^i(x — r, v) sin A n r dr 



dz/. 



(28) 
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where 



A n (x, s) = A(x, s) + — (V(s) - A > i(s, s)) sin A n (x - s) 



Note that 



0. 



3.2 Some estimates 

The estimates in the following Lemma El are used also in the Appendix, where 
we need to keep track of the dependence on T G [0, if], hence on S in the 
corresponding interval [0, H]. So, in this lemma dependence on S is explicitly 
indicated. For this, we recall that the function Y n (x) = Y n (x; S) is defined 
for < x < S < H. Furthermore we recall that also yo and y\ depend on 
S: y = yo(S) y\ = yi(S). Using the asymptotic formulas (j7|) and Gronwall 
inequality we get: 

Lemma 8: Let H > 0. There exists a number M = M= such that for every 
x G [0, S] C [0, S] t/ie following inequalities hold: 



The proof is in the Appendix. 
Using Corollary [2] we get: 

Theorem 9: Let S > n be fixed. There exists N > such that the sequence 
{Y n (x)} n >N is a Riesz sequence in L 2 (0, S). 

This suggests that we can use Bari Theorem, as in [13], in order to prove that 
{Y n (x)} n >i is a Riesz sequence. This is an intermediate step we shall need 
below, but we note that the sequence to be studied is ( TTg|) which, written in 
terms of Y n (x) gives the following sequence of functions of the variable t: 



Y n (x;S)\<M, 



Y n (x; S) — yo(S) cosnx| < 



M 



(29) 



n 




L(t) 



1 



</Q(M(s)) 



Y n (s) ds. 



The transformation from ip(t) G L 2 (0,T) to {C^){x) = ip (M(x)) in L 2 (0, S) 
(where T = M(S)) defines a bounded and boundedly invertible transforma- 
tion of L 2 (0, T) onto L 2 (0, S) so that we can equivalently study the following 
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sequence in L 2 (0, S): 

^' n (0)£c(x,s)Y n (s)ds^ (30) 

where 

C(x, s) = N(M(x) - M(s))e 7V ' (0)A/(s)/2 Q" 1/4 (M(s)) (31) 

and we must identify a value of S such that the sequence in (I3UI) is Riesz in 
L 2 (0,S). We shall prove: 

Theorem 10: The sequence / T50j) is a Riesz sequence in L 2 (0,S) for every 
S>n. 

Using ( 124)) and L'(0) = we see that the sequence ( !T8|) is a Riesz sequence 
in L 2 (0,T) for every T > M(ir), i.e. for every T > T where T solve 

° ^Q(s) &s = ix. (32) 

This is the statement in Theorem [5j 

Thanks to the fact that C(x, x) is strictly positive, we can equivalently 
study the sequence 

j&(0) jf B{x, s)Y n (s) da} , S(x, a) = . (33) 

The computational advantage of this transformation is that now B(x,x) = 1. 

3.2.1 The proof of Theorem [10] 

In this proof we work on a fixed interval [0, S] (hence in a fixed interval 
[0, T] so that we don't need to take track of the dependence of Y n (x) on S. 
We shall need asymptotic estimates for several different sequences. In order 
to simplify the notations, M will denote a number which does not depend 
on the index of the sequences but which in general depends on the interval 
[0,S] we are working with. Analogously, {M n (x)} will denote a sequence of 
functions which is bounded on [0, S]. The bound M can depend on S. These 
constants and sequences will not be the same at every occurrence, without 
any possibility of confusion. 
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We introduce 

K(s)=V(s)-A, 1 {s,s) (34) 



so that A n (x, s) in formula (|28|) is 



A n (x, s) = A(x, s) + —K(s) sin \ n {x — s) . 

Furthermore, we introduce the notation Z n (x): in ( 1331) : 

Z n (z) = n (O) / s)y n (s) ds . (35) 
Jo 

We shall proceed in parallel with the sequences {Y n (x)} and {Z n (x)} in 
order to prove the following result which, as we noted, implies Theorem 

Theorem 11: Let S > ir. Both the sequences {Y n (x)} and {Z n (x)} are Riesz 
sequences in L 2 (0, S). 

The proof require several steps. 

Step 1: linear independence. We first prove that if {Z n (x)} is linearly 
dependent then (Y^x)} is linearly dependent too; and then we prove linear 
independence of {Y n (x)}. 

If 



^a n Z n (x) = i.e. / B[ 

n=l ^0 



X, s 



K 



n=l 



ds = 



then we have 

K 
n=l 

because B(x, s) is smooth and B(x,x) = 1. 

Now we prove that {Y n (x)} is linearly independent. We prceed by con- 
tradiction: let K be the first index for which 

K 



Y,anY n {s) = 0. (36) 



n=l 



Note that K > 2 since Y x {x) ^ 0. 



Viscoelastic string under external traction 



18 



Using equality (J25]) we get 



K 



o = 5> n :C( 

n=l 

So we have also 



x 



^ \ 2 n a n Y n - / A(x,r) ^A^a n y n (r) dr 

n=l 70 \n=l / 



A" 



^A^a n F n (x) = 0. 



n=l 



This equality and (136|) contradict the definition of K since they give: 



K-l 



^T(Ai--A> n y n (x) = o. 



n=l 



Step 2: a new expression for Z n (x). We integrate by parts in (T2"g|) as 
follows (the definition of is in 



Y n (x) = g n (x) - I A(x, s)Y n (s) ds - ^- / sin \ n (x - s)K(s)Y n (s) ds 
'o * n J 



JO 



9n{X) 



d 

A u{x — s, v)—r- cos \ n s ds 



ds 



YJv) dv 



PX ^ PX 

/ A(x, s)Y n (s) ds — — / sin X n (x — s)K(s)Y n (s) ds 
Jo A n J 



+-zMJx) 



Let R(x,s) be the resolvent kernel of —A(x,s). Then we have (recall 
A n ~ n) 

\ 1 r 

Y n (x) = g n (x) t / sinA„(x - s)K(s)Y n (s) ds 

\i JO 



+ / R{x,s) 
Jo 

+—M n (x) . 



1 f s 

9n(s) - — K(v)sm\ n (s - u)Y n (u) dv 
\i Jo 



ds 



n- 



(37) 



Using smoothness of R(x, s) and R(x, x) = (a consequence of A(x, x) = 0), 
we can integrate by parts twice and we see: 
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T-> / \ f \ A M "( X ) 

R(x, s)g n {s) ds = — . 

Using Lemma [8], we see that the same holds also for the double integral 
(exchange the order of integration and integrate by parts once) so that 



YJx) 



9n{X) 



;MJx) . 



l r 

— / sinA n (x - s)K(s)Y n (s) ds 
We replace this expression in (1351) and we see that: 

f x 4>' (o) f x 

Z n (x) = 4>' n (0)y / B(x,s) cos X n s ds + yi-\ — / B(x, s) sin X n s ds 

JO *n JO 

0n.(O) 



A n J 



B(x, s) sin X n (s — v) ds 



K(u)YJu) du 



i r 

+— / B(x,s)M n (s) ds 

Jo 



(38) 



Integrating by parts, we see that on every interval [0, S] the functions 
Z n (x) are sum of a term of the order 1/n plus: 



r <t>' (o) r r 

0^(O)2/o / B(x, s) cos X n s ds = yo-\ — sinA n x- / B t 2 (x, s) sin X n s ds 
Jo x n |_ Jo 



(39) 

A further integration by parts of the last integral and the estimates (El) 
and (jHJ) show the following result: 

Lemma 12: Let S > 0. There exists M such that for every n and every 
x G [0, S] we have 



\Z n {x)\<M, 



\Z n (x) - y \ - smnx\ < — . 

V 7r n 



(40) 



Step 3: we prove that w-independence of {Y n (x)} on L 2 (0, S), S > it, implies 

that of {Z n (x)}. We proceed by contradiction. Let {a n } be a sequence such 
that the following equality holds in L 2 (0, S): 



22a n Z n (x) = 



(41) 



n=l 
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The fact that {Z n (x)} n> ^ is a Riesz sequence in L 2 (0, S) implies that {«„,} G 
I 2 . In fact, we have more: 

Lemma 13: // ( fjTp holds in L 2 (0, 5) (any S > 0) then the series 

^2<f)' n (0)a n smnx (42) 

converges in L 2 (0,S). So, If ffl\ ) holds and if S > n then there exists 
{~in\ G I 2 such that 

In 



Proof . The second statement follows from the first one, since the sequence 
{sin A n x} n> Ar is a Riesz sequence^ in L 2 (0, S) when S > n so that convergence 
of ( 1421) implies {4>' n (0)a n } G I 2 . So, we prove that equality ( flll implies 
convergence of the series (j42]) . 

We prove the first statement. Using f l38j) . equality f )4T|) can be written as 

V^^« n / / B(x,s) S m\ n (s-v)ds K{v)Y n {v)dv 
n=1 K Jo Uu J 

B{x, s)M n (s) ds-yj B(x, s) V ^-a n sin \ n s c 

t^^nJo JO Vn^t An 



+oo „ x > 

y~^0^(O)a n / -B(x, s) cos A n s ds 



(these equalities technically have to be intended as finite sums, and the equal- 
ity holds in the limit if it happens that the series converges. But, every series 
a part possibly the one on the right hand side is clearly convergent, so that 
the series on the right side has to converge too). Using we write the 
series on the right hand side as 

Vo I > 'o. n ' rn ^' sin\ n x I - y I B 2 {x,s) I } a n ^ n }^ sm\ n s ds I ds 



'+00 ^ /qn \ „ x / +00 

y^ a n-\ — sinA„x \-yo / B >2 (x, s) I V\ 

v n=l K J J \n=l 



A, 



We deduce for this that every term, a part possibly the first series above, 
can be differentiated in L 2 and that the differentiated series converge. Hence 



1 here in fact N = 1 but we don't need this. 
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also the first series can be differentiated termwise and 

' +00 



yo 



y^q ra <// n (0) cos A n 



x 



,n=l 



converges to a square integrable function, so that {0^(O)a n } G I 2 since 
{cos A n x} n> Ar is a Riesz sequence in L 2 (0, S) when S > ir. 1 

Consequently we have 



S(z,s) 



n=l 



ds 



The fact that B(x,s) is smooth and B(x,x) = 1 easily implies 

+00 

^ 7n r n ( s ) = o. 



n=l 



Arguing by contradiction, we see that if [Y n {x)} is ui- independent, the 
same holds of {Z n (x)}. So, in order to prove that {Z n (x)} is w-independent, 
which is our final goal, it is sufficient to prove ^-independence of [Y n (x)}. 
This we do in the last step. 

Step 4: the sequence {Y n (x)} is w-independent. Let the sequence {7^} 
satisfy 



+ OC 



J2lnY n (x) =0 

n=l 

in L 2 (0, S) (so that {7 n } G I 2 ). Using ([28]) we get 



(43) 



A re (z, s)y n (s) ds 



V — 



A; 



A.nfx, v)YJu) du 



+ / cos A„(x — s) 
Jo 



A n (s,s)y n (s) 



A m (s,i/)Y re (^) dv 



ds 



We prove that we can termwise compute the derivative of the series on the 
left hand side. As above, it is sufficient that we prove that the series on the 
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right side can be differentiated termwise. This is clear for the second series. 
As to the first, it can be written as 

PX PX 

V^7n / A n (x, s)Y n (s) ds = V]7n / A(x,s)Y n (s) ds 
Jo Jo 

7 [ x 

-V"t^ / K{s) sin X n (x - s)Y n (s) ds . 

*n Jo 

As to the second series we have the following result, used also below: 

Lemma 14: The series (XXTnMn) Jq K(s) sin X n (x — s)Y n (s) ds} converges 
in L 2 (0,S). 

Proof . We note that 

7 f x 

V" / K(s) sin X n (x - s)Y n (s) ds 
— A n Jo 

7 f x 

= V"-^ / K(s) \Y n (s) - ?/ cos X n s] sin X n (x - s) ds 
\% Jo 

E7n r 
— I K(s) sin X n (x — s) cos X n s ds . 
X n Jo 

The derivative of the first series converges uniformly, thanks to inequal- 
ity (12 91) . Using trigonometric formulas we see that the last series is 



I K(s) dsj^^sinX n x + ^^- K((x - r)/2) sin X n r dr 

whose derivative is L 2 -convergent. I 

We sum up: Yln=i 7™ cos ^ nX e W 1,2 (0, S) and, proceeding as in previous 
steps, we see that the following lemma holds: 

Lemma 15: Let equality (^STJ hold for the sequence {7 n }- Then, there exists 
{c n } G I 2 such that 7 n = a n /X n . 

Now we start a bootstrap argument: 

Lemma 16: The series ^2(o~ n / X n )Y n (x) is uniformly convergent to a W 1,2 (0, S) 
function, and its derivative can be computed twice termwise. 
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Proof . We insert 7 n = o n j A n in (H5|) and we use the representation 
get 



to 



71=1 n=l " 



s) sin X n (x — s)Y n (s) ds 



+00 

£ 

71=1 



On 

AH 



sin XJx — s) 



A j i 1 (s,v)Y n (v) dv 



ds 



We compute a first derivative termwise and we get 



n=l A " n=l Xn J 



s) cos\ n (x — s)Y n (s) ds 



+00 

EO~n 



n=l 



n JO 



An ( x ~ r ) v ) cos A n ?" dr 



KM dz/ 



Clearly we can compute a second derivative of the last series. We prove 
that also the first series can be again differentiated termwise. In fact, formal 
differentiation gives 



(+00 \ +00 „ x 

71=1 n / 71=1 



s) sin A^fx — slYnfs) ds . 



The first series is uniformly convergent and the second one is handles as in 
Lemma HU 

The previous computations prove L 2 -convergence of X] r t=i( <J ™/ '^n)g'n{x) 
so that, as above: 

Lemma 17: If condition §£3j) holds then there exists a sequence {5 n } G I 2 such 
that 

0~n 



In 



This lemma combined with ( 125]) shows that the second derivative of the series 
in ( )43|) can be computed termwise, and of course it is zero. So, from (j25|) we 

get 



+00 px 

S"5 n Y n {x)+ / A(x,r) 

n=l Jo 



+00 



n=l 



+00 



dr = i.e. 5 n Y n (x) = . 



71=1 
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This equality, combined with (|4~3|) . gives 



+ 00 , X 

E 1 - T2 w*) 



0. 



n=2 ^ 



This argument can be repeated till we remove N first elements from the series 
in (|43p . Using 1 — 1/A| ^ for n large we conclude that the series in f )43|) 
is in fact a finite sum. But, then all its coefficients have to be zero, since we 
proved linear independence of the sequence {Y n (x)}. 

This completes the proof of Theorem [TTJ hence also of Theorem [51 
Using this result, we can now study the identification problem. 

4 Source reconstruction 

It is now convenient to go back to the full notation z n (s; t) and we recall that 
z n (s;t) is defined for < s < t. So, the solution of Eq. ([T]) with boundary 
control f{t) =0 is given by (see ( TT4~|) ) 



W (Z> f ) = J2 / z n (t-s- t)g(s) ds, b n = (b, 4> n ) = / 6(O0n(O d^ • 



n=l 





The output y(t) 



w x (0, t) is 




(44) 



We prove that we can compute the output for every t, i.e.: 



Lemma 18: The function t — > 77 (t) is continuous. 



Proof. 



The proof is in the Appendix. 

Now we choose g(t) of the special form 




(45) 
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Here f(t) is locally square integrable and this is achieved by taking 



a(t) = f(t)+ / N'(t-s)f(s)ds 



(use N(0) = 1.) Let us fix any T >Tq (To is specified in Theorem [5]) and let 
us note that 

V( T ) = (J2Mx)<P' n (0) J o z n (T-r;T) ^ N(r - s)f(s) ds dr, b n (x) 



+oo „t 
,n=l Jo 



<j>' n (0) / N(s-r)z n (r;T)dr 



ds, b(x) 



(the crochet denotes L 2 (0, 7r)-inner product.) We recall z n (r;T) = z n {r) and 
we compare this formula with (1161) . We see that the left side of the crochet is 
nothing else then w(T) when f(t) is the boundary control. Hence, for every k 
there exists a suitable fkif) e L 2 (0,T), i.e. a suitable o~k{t) e L 2 (0,T) such 
that the corresponding output r){k){t) is such that 

V(k)(T) = ((f> k ,b) = b k , 

the k — th Fourier coefficient of £>(£) respect to the orthonormal basis {<f>k(Q} 
(or, if we whish, to any prescribed orthonormal basis, thanks to Theorem [5j) 
So, &(£) is given by 

+oo 

b(0 = Y, 7 i(k ) (T)MO- 

n=l 



Remark 19: At first sight it might seem that the previous method is similar to 
the one developed in fl&]/ . In fact, it is not the same idea and, as has to be 
expected due to the time varying coefficient P(t), it is less efficient. In fact, 
the algorithm in FTSj/ does not really use o~(t) which has to be known, smooth 
and with a(0) = 1, but fixed once and for all. Using this fixed input the output 
y(t) is measured for every t G [0,T] and then a sequence of algorithms applied 
to this simple observation produce the Fourier coefficients bk- The algorithm 
in JT^I has been extended to systems with memory, but with constant P(t), 
in but it seems that this algorithm cannot be used to study the time 
varying traction. For this reason we proposed to use the full map L 2 (0,T) 3 
a(-)^y(T). I 
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5 Appendix: The proof of Lemmas IH [H] and [TH 



In this appendix it is important to recall that every transformation in Sec- 
tion [3]Tj] is on a fixed interval [0,T] and does depend on T. Even the trans- 
formation H{t) = H(t;T) since it depends on Q(t) = P(T — t). Now we 
considere these transformation for every fixed T in an interval [0, K] so that 
it will beO<t<T<K and the pair (t, T) will belong to the triangle 

A = {(t,T) : < t < T < K}. 

We have that 



L(t) = L(t;T) = [ y/P(T - s) di 
Jo 



is bounded on the triangle A and so also is L'(t; T), uniformly for < T < K. 
As in Section 13.1. 1[ let, for every T G [0, K], 

S = L(T) = L(T; T) 

Note that S = L(T; T) so that < x = L(t; T) < S < E = L(K, K) and 
that the triangle A is transformed to 

A = {(x,S) : < x < S < E}. 

So, M(t; X) and M'(t; X) and A(x, r; X) and their derivatives are bounded 
of A. We recall that yo and y\ do depend on S. So we shall write 

9n(x) = g n (x; S) = y (S) cos X n x + — ^— y x {S) sin X n x 
and that yo(S) and yi(S) are bounded on [0, E]. 



Proof of Lemma [HI The first inequality in ([29]) follows from Gronwall 
inequality applied to fT28|) . So, we must prove 

\Y n (x;S)-g n (x;S)\ < — (46) 

n 

(the constant M depends on E.) 

In order to get inequality ( H6|) we integrate by parts once the last integral 
in (128]) and we add and subtract g n (x; S) to the first integral. Then we 
integrate by parts so to get a factor 1/A n . We obtain 
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Y n (x) - g n (x; S) = 

- [ X A n (x,r;S)[Y n (r;S)-g n (r;S)] ds 
Jo 

A n (x, r; S)g n (s; S) dr + —M n (x; S) 

K 

where M n (x; S) is a function that we don't need to write down explicitly, but 
such that \M n (x, S)\ < M for every n and (x, S) G A. 

The require inequality follows from here, integrating by parts the last 
integral and using Gronwall inequality, i 

The proofs of Lemmas|4]and [TBI The idea of the proofs of these lemmas 
is similar. In order to prove Lemma H] we have to consider the first series 
in (fl4|) while for Lemma [18] we have to consider the series in ( 144)) . The 
common feature is an integral of the general form 

/ V(t,r)z n (r;T) dr 
Jo 

L(t;T) 

V(t; M(s- 5))e 7V (0)m(s;5) P- 1/4 (M(s; S) - T%{M{s; S); T) ds . 
Now we use 

L(M(x; S);T)=x, Y n (s; S) = C„(M( S ; S); T) 
and finally we get an integral of the form 
"an" era scritto "and" 




/ V 1 (x,s;X)Y n (s;X)ds. 
Jo 

The integral in Lemma H] is obtained when V(t, r;T) = N(t — r) and then 
Vi(x,s\X) = C(x,s;X). The integral in (jHJ) is obtained with V(t,r) = 
g(t — r). Note that C(x,x;X) ^ while g(0) = and so in the proof of 
Lemma [TS] we shall have Vi(s, s; X) = 0. 
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We present now the proof of Lemma H] and we leave the similar proof 
(based on inequality |4"0"|) of Lemma [TS] to the reader. 

We transform the variable t G [0, T] C [0, K] in the series f fl4|) to the 
variable x = M(t) G [0, S] C [0, S] as shown above and we consider the 
L 2 (0, 7r) norm of the resulting series, for each S G [0, S]. Using the fact that 
{</>„(£)} is orthonormal in L 2 (0, 7r) and definition f[3"Tj) . the square of the norm 
can be written as 



E 



n=l 



K(x,v;S)\4>' n (0) I C(v,r;S)Y„(r;S)Ar\ iv 



E 



n=l 



A(x, ^; S 1 ) (C(^, z/; S)Z n (v; S)) du 
A(x, !/; S) = (f (M(x; S) - M(u; S)) M\v- S)) 
Using fHUj) we see that it is sufficient that we note 



E 

71=1 



A(x, v\ S)C(u, v\ S) sinns 



< 



\K{x,v-S)C{v,v-S)\ 2 ds. 



The proof is now finished since the right hand side is bounded for < 
s < S < E. i 
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